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We introduce the stochastic domino cellular automaton model exhibiting avalanches. Depending 
of the choice of the parameters, the model covers wide range of properties: various types of expo- 
nential and long tail (up to inverse-power) distributions of avalanches are observed. The stationary 
state of automaton is described by a set of nonlinear discrete equations derived in an exact way 
from elementary combinatorial arguments. These equations allow to derive formulas explaining both 
various exponential and inverse power distributions relating them to values of the parameters. The 
exact relations between the state variable of the model (moments) are derived in two ways: from 
direct arguments and from the set of equations. Excellent agreement of the obtained analytical 
results with numerical simulations is observed. 
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I. INTRODUCTION 

There is a significant interest in constructing simple 
stochastic models with avalanches reflecting properties 
of various natural phenomena, but only part of them 
have the advantage of being solvable. One way to make 
models analytically tractable is to use the mean-field ap- 
proximation. However, requirements of exact results can 
force further simplifications in a model. For example, 
in terms of self-organized criticality, in spite of the fact 
that the mean field theory of such critical phenomena 
was already proposed [lj], Dhar made further simplifica- 
tions for distinguished BTW model @ (a unified view 
of mean-field picture of stochastic self-organized critical 
models is presented in Q). We point out also that con- 
sidering structure of abelian algebras leads to analytical 
results in some sand-pile and related stochastic models 
d, d| . Solvable simple models can also be constructed in 
the field of directed percolation; for an extremely simple 
one, see @. 

Another way towards solvability is to consider stochas- 
tic properties as non-essential and study deterministic 
models; see such an approach for sand-pile in [7j. In the 
context of deterministic cellular automata, an analytic 
approach related to solvability is investigated as part of 
the theory of integrable systems and may employ sophis- 
ticated methods. |cj-[l3| . We underline that construction 
of integrable cellular automata is definitely not an easy 
task. 

Our aim here is to follow the ideas of analytic approach 
and apply them to construction of stochastic automata. 
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In the article we propose and analyse in a direct elemen- 
tary way the Random Domino Automaton (RDA) - a new 
slowly driven systems exhibiting avalanche phenomena. 
We prefer here an elementary self-contained approach to 
the description of the automaton and we make no use 
of applicable Markov processes terminology. We stress 
here, all results, except of equation (|28l) (see [l4| for de- 
tails), have elementary derivation in the text below. On 
the other hand, in spite of its simple formulation and be- 
ing analytically tractable, automaton covers wide range 
of behaviours depending on the choice of the parame- 
ters. The application of the RDA model for studying Ito 
equation is investigated in our parallel papers [15|, [l6( . 

An inspiration for defining the rule for the domino au- 
tomaton comes from very simplified view of earthquakes. 
It corresponds to two tectonic plates moving with rel- 
ative constant velocity. The wedge may be irregularly 
rough, and relative motion can be locked in some places, 
producing stress accumulation. Beyond some threshold 
of stress, a relaxation took place. The size of relaxation 
depends on the nearby accumulated stress. RDA is in- 
spired by earthquakes; however, a direct reproducing of 
realistic-like behaviour of such complicated phenomena 
is obviously beyond its scope in the present form. Never- 
theless, construction and analysis of models is one of pri- 
mary aims in geosciences [l7l [l8j . There are many very 
simplified cellular automata models focusing on specific 
features of the investigated behaviour in the field. Here 
we point out a sequence of papers fl9rl2lj ]. where some 
interesting cellular automata models were presented. 

Finally, as an unexpected property we mention an in- 
teresting link between stochastic cellular automata and 
integer sequences (see [22j or The On-Line Encyclopedia 
of Integer Sequences). It is known, how to obtain Cata- 
lan numbers out of the bond directed percolation on a 
square lattice |23[. Our cellular automaton in specific 
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case leads to Motzkin numbers (for details see [HI). 

Plan of the article is as follows. In Section [TT| we in- 
troduce a definition of random domino automaton with 
rebound parameters. Section IIIII contains full combina- 
torial derivation of equations for the distribution of clus- 
ters, which describes stationary state of the automaton. 
Several exact formulas - balance relations, average size of 
clusters and avalanche, general formula for all moments 
- are displayed. We finish with introducing special form 
of rebound parameters, which leads to two distinguished 
special cases studied in details in next Section HVl These 
cases correspond to exponential and inverse-power distri- 
butions as shown below. The obtained analytical results 
are compared with simulations. Section |V] gives a resume 
of results presented in the article as well as directions for 
future work. Appendix [X] contains derivation of comple- 
mentary set of equations, describing empty clusters of 
the automaton. 



II. DEFINITION OF RANDOM DOMINO 
AUTOMATON 

In the random domino automaton model, the space 
consists of discrete number ./V of cells on a line, and we 
assume periodic boundary conditions. Each cell may be 
in one of two states: empty state, when it is empty, or oc- 
cupied state, when contains a ball. The evolution of the 
automaton is given by the following update procedure 
performed in each discrete time step. A ball is added 
to the system to the randomly chosen cell and we as- 
sume each cell to be equally possible. If the chosen cell 
is empty, there are two possibilities: it becomes occupied 
with probability v or the ball is rebounded with probabil- 
ity (1 — v) leaving the state of the automaton unchanged. 
If the chosen place is already occupied, there are also 
two possibilities: the ball is rebounded with probability 
(1 — [i) or with probability p, the incoming ball triggers a 
relaxation. By relaxation we mean: balls from the chosen 
cell and from all its adjacent occupied cells are removed. 
Thus, the relaxation produces an avalanche of size equal 
to the number of cells changing their state. Then the 
update procedure repeats in the next time step. All pos- 
sibilities are shown schematically on the diagram below. 

• 
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An avalanche is represented by a symbol •. An example 
of relaxation of the size three is presented in the diagram 
below. 
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The name 'domino automaton' comes from the fol- 
lowing interpretation. Occupied cells are represented by 
standing domino blocks and empty cells are represented 
by extra space between them. If the incoming ball hits an 
empty space, a domino block is added there with prob- 
ability v. If the incoming ball strikes a domino block, 
with probability p it falls down the chosen block and all 
its adjacent neighbours (on both sides) up to the gap 
(empty cell). Then thefallen dominoes are removed and 
the procedure repeats. 

Presented defining rules of the domino automaton can 
be easily modified leading to various extensions of the 
system. We mention some possibilities (like geometry, 
capacity of cells, many kinds of balls, different triggering 
rules) in section [Vj However, the aim of this work is to 
provide detailed description of the "core" case described 
above. 



III. EQUATIONS OF RANDOM DOMINO 
AUTOMATON 

A. Notation 

A state of the automaton is defined by altered se- 
quences of occupied and empty cells. A sequence of i 
consecutive occupied cells is called the cluster of length i 
(shortly i-cluster); a sequence of i consecutive empty cells 
is called empty cluster of length i, where i = 1, 2, . . . , N. 
The fixed size of the lattice N is assumed to be finite, but 
big enough to make limit for the size of clusters negligi- 
ble. Denote the number of i-clusters by rij, the number 
of empty i-clusters by n®, the total number of clusters by 
n, and the density (a number of occupied cells divided 
by the number of all cells N) by p. Then, it follows that 

n=1 ^2 n i = ^2 n i and p = -^^2 n > L C 1 ) 

i>l i>l i>l 

For higher densities p, the probability of relaxation is 
also relatively higher, and the density is more likely to 
decrease. For smaller densities, triggering of an avalanche 
is relatively less probable and the density tends to grow 
up. Hence, the variations of the density (see FigQJ are 
subjected to 'v-shape' potential. The behaviour of the 
automaton is described below under the assumption that 
it is in a quasi-equilibrium state and the variables used in 
the equations, like density and others, are average values 
and do not depend on time. The applicable description of 
the system as a Markov process is postponed to another 
paper. 
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FIG. 1: Examples of simulation results for time series of den- 
sity p(i) of the ID domino automaton in case p/v = 1 with 
lattice size N = 500 and in case fi/v = with lattice size 
N — 4000. The parameter i numbers avalanches. 



We consider the following settings for rebound param- 
eters. The coefficient p = pi does not depend on the 
position of the chosen cell, but may depend on the size i 
of the chosen cluster. The parameter v is fixed to be con- 
stant. There are two cases of special interest (discussed 
in detail in Section HV)). One with p — (3 — const refers 
to equal probability of provoking an avalanche for each 
occupied cell, while case p, = /i,- = j, where 6 = const, 
refers to equal probability of triggering an avalanche for 
each cluster, irrelevantly to its length. These two respec- 
tive cases are related to exponential and long-tail distri- 
butions of clusters. 



B. The balance equation for p. 

The stationarity condition requires the flow-in to be 
equal to flow-out, hence there must be a balance between 
losses and gains in the number of occupied cells. In a 
single time step, the number of occupied cells may stay 
unchanged (when ball is reflected), or may increase by 
one (an empty cell become occupied), or decrease by i 
(in a case of an avalanche of size i). The expected value 
of increase of density is equal to the probability that an 
empty space becomes occupied, namely v(\ — p). The 
probability of relaxation of a size i is /ij^. Since any 
possible size i can trigger an avalanche, the stationarity 
gives 



N 



i>l 



v 



C. The balance equation for n. 

In a single time step, the number of clusters may in- 
crease only when a new 1-cluster is created. The chance 
is proportional to the number of interior cells in empty 
clusters having the length three and bigger, as depicted 
in the diagram below. 

(i— 2) — interior 



Hence, the probability is 
V- (i-2)n?i 



j>3 



N 



(3) 



Losses in the total number of clusters may appear in 
two ways: by joining of two clusters (separated by an 
empty 1-cluster) and by triggering an avalanche. In the 

first case, the probability is ~ c-J-, in the second it is ~ 

Si>i P-i^w- Hence, the balance equation for the number 
of clusters is of the form 



(l-p)N-2n = J2- 



(4) 



i>l 



D. Balance equations for ms. 

We underline the assumption used in the derivation 
of the equations below: clusters are distributed indepen- 
dently, by which we mean that the length of the "next" 
cluster does not depend on the length of the " previous" 
one. In other words, our investigations are done up to 
the order of clusters (and order of empty clusters). The 
presented approach may be regarded as a generalization 
of percolation, where subsequent cells are treated as in- 
dependent. Below, in subsection IIV Al we compare both 
approaches numerically. 

To write down equations for the numbers of clusters 
of length i, i.e., for ma, we consider all possibilities of 
losses of such clusters as well as creation of them, and 
next we claim that on the average the gains and losses of 
respective values compensate each other, as required by 
the stationarity conditions. 

Losses. There are two ways to destroy an «-cluster: 
by enlarging and by provoking the avalanche depending 
on the cell where an incoming particle is thrown. 



(a) Enlarging. For any cluster there are two cells adja- 
cent to its ends, so the probability is 



(2) 



2v- 



N' 
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If the single empty cell is between two clusters of the 
length i it is counted twice - it decrease the number of 
i-clusters by two. 

(b) Relaxation. In this case it is enough to knock out 
any of the occupied cell of the cluster, so the probability 



is 



f'i 



N 



Gains. There are in general two possibilities to create 
z-cluster: enlarging (i — l)cluster and joining two smaller 
clusters. 

(a) Enlarging. Case i — 1 was already considered and the 
probability is given by formula For i > 2 enlarging 
a (i — l)-cluster to the size i is possible if the adjacent 
empty cluster is of a size bigger then one. Hence the 
probability is 



2v- 



i-iE 



i>2 



N 



2v 



N 



where the multiplier 2 counts left and right cases, 
(b) Joining two clusters. Two smaller clusters: one of 
size k G {1, 2, . . . , (t — 2)} and the other of size (i — l) — k 
will be joined if the ball fills an empty cell between them. 




The probability is proportional to the number of empty 
1-clusters between A:-cluster and (i — 1 — fc)-cluster, hence 



o »~ 2 
n" \ - n k rii 



N f—* n 

k=l 



1-fc 



The dot in the multiplication above underlines the inde- 
pendence assumption for the order of clusters. The last 
formula introduces also an extra quadratic nonlinearity 
into the system; so far, nonlinearity was present through 
n (and also n\). 

Finally the following set of equations for n^s is derived 



1 



ni 



n- 2 = 



Ei +2 

V 

2 

2^ + 2 

V 

1 



((1- p)N -2n + nf) 



1 - 



X 2r^i l-£ + n?£ 



(5) 
(6) 

(7) 



fc=i 



for i > 3, where n = J2i>i n i an d P = j? J2i>i in i- 

Variable n\ can be cancelled from the above set by 
considering the balance for empty clusters and deriving, 



in analogous way, the respective set of equations for vari- 
ables n° (see Appendix [A] for details)). The result is 
given by equation (|A1[) . namely 



2n 



(8) 



The above, substituted into equations ©-(JT]), form a 
closed set for variables {n\,n-2, • • •}• Then, equations 
(jAlj) and ([A2jl allow us to find n° k for all k = 1, 2, . . .. 

In conclusion, we obtained closed, nonlinear set of 
equations for the distribution of the clusters in the cellu- 
lar automaton. 



E. Moments 

The balance equations (|2|) and ((U (as well as equations 
for higher weighted moments of n,i) can be also obtained 
from the above set ©-jZ]). 

For m, i = 1,2,..., we define a moment of order 7 by 



at E'' 



i>l 



(9) 



Then, the density p is equal to the first moment mi, and 
the normalized number of clusters is equal to the zero 
moment tuq. We define also weighted moments m 7 as 
follows 



~ N ^ v 

i>l 



(10) 



In this notation, equation © takes the following form 
_ 2m N 



3 2 mi 

m 



(11) 



Removing denominators from the set of equations (j5])-(jTJ) , 
multiplying both sides by i z and performing summation 
with respect to i, one obtains 



E 

i>l 



M« -2 + 1 



2^2i z rii = (1- p)N -2n + n?+ (12) 



;>i 



2(1 - ^) + X )* n ' + S EE 



i>l 



i>3 fe=l 



Finally, after using formula (|11[) , changing the order of 
summation (in indices i and fc), introducing variable j by 
formula i = k+l+j, expanding binomials and performing 
sums, one obtains 



m z+1 = 1 - mi - 2m + 2 E ^ ■ 1 



(13) 



2 ^7 « v+fV 

3m + 2mi ^ V + Z / ' TOp ' 
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The above equation for z = and 2 = 1: 

rh\ = 1 — mi — 2mo, 
m 2 = I — mi, 



(14) 
(15) 



are the balance of n equation ((4]) and the balance of p 
equation ^ , respectively. Equations of moments for sub- 
sequent values of z allow to choose particular forms of 
rebound parameters and isolate interesting cases as pre- 
sented below for formula (|18l) . 

Definition of moments ^ leads to particularly neat 
form of expressions for the average size of a cluster 



< i >= 



Si>l 



and the average size of an avalanche 



m : 



< W > = 



m 2 
fhi 



mi 



J2 l >il J -i n i i ™-i l-mi-2m 
For rebound parameters in the form 
S 



Mi = — . 



const., 



(16) 



(17) 



(18) 



where 5, a and 9 = £ are constants, equations (|14p and 
(fT5")) take the form 



mi + 9m 2 -a = 1, 
2mo + mi + #mi_ CT = 1. 



(19) 
(20) 



Hence, cases <r = and cr = 1, considered in details 
below, are very suitable for exact analysis. There is also 
an interesting behaviour when a > 2, but we consider it 
in another paper. 



IV. SPECIAL CASES 

A. Case p, = /3 = const: equal probability of 
triggering relaxation for each occupied cell. 

To get the value of density from equation ^ in the 
just considered case, extra relations are required. As the 
first step, consider the percolation approximation [2~i| 



c(l - p)V 



(21) 



which means that cells are treated as independent. From 
equation (JXJ) it follows that c = N. Then equation @ 

gives p = 2/(3| + ^9(£) 2 +4(l-£)). Any value of 
density p £ (0, 1) can be obtained for suitable - € (0, oo). 
The percolation approximation in the case § = 1 gives 
the value of p = ^ , which differs by several percent from 
the value < p >~ 0.3075 obtained from numerical sim- 
ulation of the automaton. This divergence appears be- 
cause the automaton rule states that there is a " coupling" 



n(i) 

5U u'* 



FIG. 2: Values of nt in case p = /? = const from simulation 
(dots) compared with values computed from equations (solid 
line) and the approximation ke n (dashed line). Lattice size 
N = 500 and p/v = 1. 



TABLE I: The average density < p >, the average size of a 
cluster < i >, and the average size of an avalanche < w > 
from simulation results, equations of the model and percola- 
tion approximation for case p/v — 1. 





Simulation 


Equations 


Percolation 


<p> 


0.3075 


0.308 


1/3 


< i > 


1.5974 


1.597 


3/2 


< w > 


2.2487 


2.252 
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"Based on the density balance. 

between the adjacent cells, since the relaxations takes out 
the whole cluster, so treating cells as independent, like 
in formula (|21j) . gives only very rough estimation of the 
density p. 

In the general case (not percolation approximation), 
fj, = ft = const, the balance equation for n gives the 
following formula for an average size of the cluster 



< i >= 



Np 



2p 



(22) 



n l-p(l + £) 
and the balance of p gives an average size of the avalanche 



< w >- 



(23) 



A comparison of values of the average density < p >, 
the average size of cluster < i > and the average size 
of the avalanche < w > in exemplary case p:/v — 1 ob- 
tained from numerical solution of equations, from perco- 
lation approximation results (based on the density bal- 
ance) , and from simulation is presented in Table HI The 
lattice size is fixed as N — 500. 

Approximate solution for rt^s in the case p, = /3 = 
const can be obtained by substituting in equation ([7]) 
the following formula 

m = fce~ 7i for i = 3, 4, . . . (24) 

where k and 7 are some constants. For i » 1 this gives 



= _ 7 _ v kn\ 

j3 n 2 ' ' \/3 n 2 
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FIG. 3: Exact values of m obtained from the model (line) 
in case fit — 8/i compared with exemplary simulation data 
(dots). Lattice size N = 4000 and 8/v = 1/4. 



TABLE II: The average density < p > and the average size 
of a cluster/avalanche < i > — < w > from simulation results 
and equations of the model for the case [n = 4, for 6 = 1/4. 





Simulation 


Equations 


< p > 


0.7926 


0.8 


< i >=< w > 


8.967 


9 



Table |TT] contains values of average density and average 
cluster / avalanche size obtained from equations and sim- 
ulation in an exemplary case where 9 = 0.25. The size of 
lattice is set to N = 4000 in order to avoid restrictions for 
size of avalanches, which are relatively bigger comparing 
to the previous case. 

A simple form of n\ = 2n/(3 + 29) together with bal- 
ance equations (|26|) and (f27|) alows to reduce the set of 
equations to the solvable recurrence. The formu- 

las for riiS are rational functions of 9 and N only. The 
explicit form of solutions and relation to the Motzkin 
numbers is presented in (T3 |. 

Figure [3] presents exact values of njS compared in log- 
log scale with the simulation result in the exemplary case 
described above. The slope of the distribution of rii ap- 
proximates inverse-power distribution in the middle part 
of a range of sizes i; for bigger values of i, values of ?Xj 
decrease quicker. This property agrees with convergence 
of moments of we use above. For value of 9 closer to 
the straight part of the plot extends for bigger values of 
i. In the limit case 9 eN _!^ nst distribution is given 

by the following power law [l4| 



The value of constant k can be found from equation |TJ . 
The above result indicates a close relation to the perco- 
lation dependence for m (see equation (l21~j) ). Approxi- 
mation (|25j) works well even for ri2, as may be seen in 
Figure [2] 



B. Case p = pi = 8/i: equal probability of triggering 
relaxation for each cluster. 

For Hi = |, where <5 = const, the balance for p - equa- 
tion © - takes the form 



1 



P = 



1 



(26) 



where 9 = —, which relates the density with the ratio of 
coefficients S and v only. Since 9 € (0,oo), any density 
p G (0,1) may be realized. The balance for n - equation 
(J4j) - is reduced to (1— p)N = (2+9)n. Together equations 
© and gl) give 



= N 



(9 + l){6- 



(27) 



The maximal number of clusters is (l+v^) -2 ^ « 0.177V, 
as obtained for 9 = \/2. 

The average size of a cluster is the same as the average 
size of an avalanche, because any cluster, regardless of 
its size, has the same chance to become an avalanche. 
Hence, both balance equations lead to the same result 

2 

< i >=< w >= 1 + -. 



1 

.3 

?.2 



(28) 



V. CONCLUSIONS. 

We proposed and studied properties of simple stochas- 
tic cellular automaton with avalanches inspired by an 
extremely simplified model of earthquakes. The main 
goal was to construct a simple model with "transpar- 
ent" mathematical structure, which would enable ob- 
taining exact results yet covering a wide range of be- 
haviours. The rule of the automaton allows to derive 
from the first principles (using elementary combinatorics) 
the set of equations describing the average values of the 
model parameters. These equations leads to equations 
for moments and to neat formulas for average cluster 
and avalanche sizes in terms of the zero and the first 
moments. 

After analysis of moment equations, we considered in 
detail two cases which differ by the form of parameter 
/i responsible for triggering of avalanches. We obtained 
exponential type and inverse-power type distributions of 
clusters within uniform framework. (The distribution of 
avalanches Wi is easily obtained from the formula u>i = 
Y^^^hi-n-i ~ Pi n ii-) The quasi-equilibrium assumption 
seems to be justified by simulation results. 

The exponential case p — j3 = const is not fully analyt- 
ically solvable; nevertheless, the closed set of equations 
leads to several exact relations. Moreover, an approxi- 
mate formula for n.; was derived, and intimate relation 
to percolation was pointed out. Hence, the proposed ap- 
proach can be regarded as an extension of percolation 
approximation results, and can provide basis for exact 
treatment of other models. The automaton rule clearly 
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states that there is a "coupling" between the adjacent 
cells, since the relaxations takes out the whole cluster, so 
treating cluster - not cells - as independent is a substan- 
tial improvement in the presented case. 

The case fa = 4 is fully solvable, and all variables 
are expressed as rational functions of ratio £ of rebound 
parameters only. The shape of the distribution of clus- 
ters rii approximate an inverse-power distribution in the 
various ranges of variable i, depending on chosen param- 
eters. In the limit case, distribution tends to ~ n~2. 
The critical density is equal to 1. The limit case of the 
presented automaton reduces the set of equations to the 
recurrence, which leads to known integer sequence - the 
Motzkin numbers (see [l4{). This result establishes a 
new, remarkable link between the combinatorial object 
and the stochastic cellular automaton. 

The first application of RDA - related to Ito equa- 
tion - is already studied in our parallel papers [HI, [l6[ . 
The model serves as a fully controlled stochastic "phe- 
nomenon" and an applicability of the reconstruction of 
the Ito equation from generated time series was tested. 
Due to its simplicity, it allows to derive exactly the suit- 
able Ito equation, and analytical results were compared 
with histogram method. The obtained results are part of 
broader studies of the privilege concept and its role for 
appearance of inverse-power distributions [26j . 

The model posses also some properties which may be 
used in seeking of applications to natural phenomena and 
checking their adequacy. For example, in the context 
of earthquakes, the discrete nature of the model makes 
possible an adjustment of positions and sizes of cells to 
the geological structure of a fault. Bending and shrinking 
of the TD grid have no influence on equations; the only 
important feature is an order of cells. In this context, 
relations describing average size of cluster < i > and 
avalanche < w > may be interpreted as relations giving 
their dependence on the density of energy on the wedge. 

The model establishes one-to-one correspondence be- 
tween distribution of avalanches (i.e. observable quan- 
tity) and respective parameter responsible for trigger- 
ing avalanches, as it was illustrated in two cases consid- 
ered above. We emphasize that the relation works in both 
directions. Fixing rebound parameters leads to various 
distributions of avalanches, but also, for any distribution 
m, one can easily find unique values of rebound parame- 
ters by simple use of the set of equations ©-(0. Thus, 
within the proposed model, one can infer about trigger- 
ing properties from "observed" statistics of avalanches. 
For an earthquake interpretation, it would give an in- 
sight into properties of microscopic mechanism of releas- 
ing energy on a fault on the basis of already collected 
data. 

Last but not least, we propose few generalizations pos- 
sibly leading towards more realistic extensions of the 
model. Each element of the presented automaton - in- 
cluding the incidence rule - can be subjected to various 
modifications. To be more specific, one can consider dif- 
ferent geometry of the array (for example, a tree shaped 



like Bethe lattice or any in bigger dimension), different 
capacities of cells and different kinds of balls distributed 
to the system. Also there are many kinds of dependence 
of energy release threshold on other parameters (on space 
position, on states of cells in a neighbourhood etc.). We 
leave these topics for further investigations. 



Appendix A: Derivation of balance equations for 
empty clusters. 

In analogy to the set of equations ©-(0) for clusters, 
one can deduce the following set of equations for the 
empty clusters. 

Losses. There are two following possibilities. 

(a) Occupation of any cell belonging to the empty k- 
cluster. The probability is just 

/fcn° 
~ v — 
N 

(b) Provoking an avalanche. An empty cluster of the 
length k is lost when one of the two adjacent clusters is 
knocked out and forms an avalanche. 



The probability is 



L k \ ~* 



in, 



It comes from the probability of knocking out of i-cluster 
times the probability that on the end of the cluster is 
empty fc-cluster, times two ends, summed for all possible 
values of i. 

Gains. Again, there are two possibilities, 
(a) Shortening an empty cluster. An empty cluster of a 
length k can be obtained from a bigger one of the length 
i > k + 1 when an incoming ball hits its (k + l)th empty 
cell 



or symmetrically, its (i — k)th cell. Hence the relevant 
gain term is of the form 



2, V H 

N 



i>k+l 



(b) Provoking an avalanche. If an incoming ball knock 
out the cluster of the length j e {l,2,...,fc — 2} sur- 
rounded by two empty clusters of the respective lengths 
I € {1,2,..., k- 1-j} and (k-j-l), 



k-i-j 
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then the empty cluster of the length k is created. There- 
fore, for k > 3 this leads to the following term 

k-2k-l-j . o o 
~ Z^ Z^ f*3 N n n 



j=l 1=1 



This is the probability of knocking out of the j-cluster 
times probability the Z-empty cluster is on the left times 
probability the empty {k — I — j)-cluster is on the right, 
for all possible values of I. 

Balance equations. The set of equations for the av- 
erage values in the quasi equilibrium is of the form 



Balance of n out of empty clusters. The set of 

equation for empty clusters can be treated by methods 
we use in Subsection MI El in order to derive relations for 
respective moments. 



Below we consider zero moment for 



which leads 



also to equation This is implied by the fact that 

the number of empty clusters is equal to the number of 
clusters (see equation (TT])). To perform sum of (|A1[) and 
I2j) for all values of k notice the following identities: 



k>l 



i>l 



>>1 



nl\k 



E 



Mi 



n '■ — ' v 



2 nl fc = l,2,(Al) 



i>k+l 



i>l 



i>k+l 

k-2k-l-j . o 

E E *>3- (A2) 



j=i i=i 

Equation (jAll) may be written as 



n (k + 2) + ± V t>1 i±im 

nl = 1 2. 5 ±p£ 1 ^ fo r k = 1, 2, 



n 



2(n-Eti 
and equation (|A2[) . valid for k > 3 is 
o (fc+j) + |E^il^ 

Thu — n ' 

l \ n - Zi=l n i> + Z/j'=l L,t=l VTf? rfi 

Thus, on right hand sides there are terms n® with j < k 
only. 



E E *? = £(i-iK = (i-p)Jv 



k>l i>k+l 



k>l 



and 



Z^Z^ Z^ ^ jv n° n° Z^ v 

fc>3i=i i=i j>i 
In the last identity we assumed that it is allowed to 
change the order of summation with respect to indices 
j and k. Thus, a sum of (|A1|) and (|A2[) for all fc leads 
directly to 



-riii. 



It is the balance equation Q already obtained. 
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